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ABSTRACT 

The  structure  of  bare-tropically  unstable  disturbances 

in  the  tropics  is  studied  with  a  two-level  quasi-geostrophic 

model.   An  Ekman  layer  is  attached  to  the  lower  boundary. 

The  equations  are  linearized  and  the  most  unstable  mode  is 

found  numerically  by  use  of  the  initial  value  technique. 

CISK  heating,  horizontal  friction  and  shear  are  introduced 

to  give  more  realistic  solutions.   Computations  are  made  for 

two  shear  zone  wind  profiles.   The  first  is  represented  by 

U",  =  U_  =  -U  tanh  y/y   and  the  second  is  a  mirror  reflection 
1    3     o        o 

of  the  first  and  is  given  by  U_  =  -U-,  .   The  model  is  run 
for  142  days  after  which  growth  rates  are  determined  as  a 
function  of  the  non-dimensional  wave  number  ky  .   The 
associated  phase  angles  and  amplitudes  for  the  wave  numbers 
with  maximum  growth  rates  are  given.   The  solutions  with 
heating  do  not  have  a  strongly  preferred  mode  of  maximum 
growth. 
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I.   INTRODUCTION 

Synoptic  scale  disturbances  in  the  tropical  atmosphere 
have  been  recognized  for  a  long  time  (Palmer  1951,  Riehl 
195*0.   This  paper  is  concerned  with  wave  disturbances  that 
occur  in  the  lower  troposphere  and  makes  no  distinction 
between  easterly  waves  and  equatorial  waves.   Starr  and 
Wallace  (1964)  showed  that  in  the  lower  troposphere  these 
disturbances  have  a  "cold  core"  structure;  the  rising  air  is 
cooler  than  the  air  which  is  subsiding.   Warm  core  structures 
at  higher  levels  with  significant  correlations  between 
temperature  and  vertical  velocity  have  been  found  by  Chang 
et  al.  (1970)  and  Nitta  (1970).   Vertical  velocities  associ- 
ated with  these  disturbances  in  general  are  rising  air  to 
the  east  of  the  wave  axis  and  subsidence  to  the  west  (Yanai 
and  Nitta,  1967). 

The  energy  source  for  these  waves  is  not  known. 
Baroclinic  instability  is  unlikely  because  of  the  smallness 
of  the  Coriolis  parameter.   Charney  (1963)  performed  a  scale 
analysis  which  indicated  that,  in  the  absence  of  condensation, 
large-scale  motions  in  the  tropics  are  governed  by  the  baro- 
tropic  vorticity  equation.   Charney  (1963)  suggested  that 
easterly  waves  arise  from  barotropic  instability.   The 
spectral  analysis  of  Wallace  and  Chang  (1969)  indicate  a 
conversion  of  kinetic  energy  from  the  mean  flow  to  the 
disturbance. 


A  third  possible  energy  source  for  these  disturbances 
is  latent  heat  of  condensation. 

This  mechanism  is  consistent  with  observations  of  warm 
core  structures  at  upper  levels  (Chang  et  al.  1970).   This 
heating  has  been  parameterized  by  Ooyama  (196*0,  Charney 
and  Eliassen  (1964),  and  Kuo  (1965)  in  such  a  manner  that 
the  condensation  heating  is  proportional  to  the  convergence 
of  moisture  in  the  Ekman  layer.   This  leads  to  a  simple 
instability  which  is  called  Conditional  Instability  of  the 
Second  Kind  (CISK).   Yamasaki  (1969)  studied  the  influence 
of  CISK  on  linearized  disturbances  in  an  atmosphere  with 
vertical  shear.   The  solutions  displayed  many  features  that 
are  observed  in  actual  tropical  wave  disturbances.   Holton 
(1971)  has  obtained  similar  results  with  a  model  that  con- 
tains a  moving  heat  source  and  vertical  wind  shear.   These 
studies  strongly  suggest  that  the  tropical  disturbances  are 
driven  by  condensation  heating.   It  is  known  that  the  CISK 
mechanism  must  have  a  pre-existing  disturbance  to  organize 
the  convective  elements  before  CISK  can  operate.   This  pre- 
existing disturbance  could  arise  from  barotropic  instability 
as  was  demonstrated  by  Bates  (1970).   Kuo  (19*19)  examined  the 
stability  characteristics  of  a  barotropic  zonal  current  and 
found  that  instability  is  possible  only  when  the  absolute 
vorticity  gradient  reverses  sign  in  the  region. 

The  purpose  of  this  thesis  is  to  extend  the  work  done 
by  Williams  et .  al.  (1971)  to  include  cases  with  CISK  heating, 
horizontal  friction  and  vertical  shear.   Since  the  Williams 


et  al.  (197D  paper  is  referenced  many  times,  it  will 
hereafter  be  referred  to  as  WSN. 


II.   THE  FORECAST  EQUATIONS 

The  quasi-geostrophic  prediction  equations  will  be  used 
as  the  forecast  equations  following  WSN.   It  is  recognized 
that  the  quasi-geostrophic  approximation  is  not  very  accurate 
near  the  Equator 3   but  Matsuno  (1966)  has  shown  with  a  baro- 
tropic  model  that  quasi-geostrophic  motions  are  identifiable 
near  the  Equator  if  the  wavelength  is  not  too  long. 
Yamasaki  (1969)  has  found  that  the  error  is  small  at  20 
degrees  latitude  in  a  study  with  a  baroclinic  model. 

The  quasi-geostrophic  equations  are  simplified  through 
the  introduction  of  a  two-level  model  which  was  also  used 
by  Bates  (1970).   The  equations  are  linearized  and  solved 
numerically  by  the  initial  value  technique.   The  growth  rates 
and  wave  structures  are  then  computed. 

A  simple  two-level  model  is  constructed  by  dividing  the 
atmosphere  into  four  layers  with  a  constant  pressure  differ- 
ential of  Ap/2  (Fig.  l)s  numbered  0  to  U    from  top  to  bottom. 
The  boundary  condition  at  p  =  p,  is  u  =0,  which  prohibits 
vertical  propagation  of  energy.   Charney  (1969)  has  shown 
that  such  energy  propagation  is  unlikely.   At  p  =  p  the 
boundary  condition  is  u  =  to  where  w  is  the  frictionally 
induced  u  at  the  top  of  the  Ekman  layer.   Charney  and 
Eliassen  (19*19)  derived  the  following  expression  for  u 


A   1 
we  =  -p^g  (gf-)   sin  2c^lj   »  t2*1) 


where  A   is  the  kinematic  eddy  viscosity,  a  is  the  surface 
inflow  angle,  and  Cn  is  the  surface  geostrophic  vorticity. 


p  Level 

p   0 

h ! 

p.  +  Ap   1 2 


t 


h. 


3 

U)=C0 

P0  9 1 


Fig.  1  Two-level  prediction  model 

It  has  been  shown  by  Holton  et  al .  (1971)  and  Chang  (1973) 
that  this  expression  is  reasonably  valid  as  far  equatorward 
as  10  degrees,  but  may  be  invalid  when  close  to  the  Equator 
or  south  of  the  ITCZ.   The  surface  geostrophic  vorticity  can 
be  approximated  by 

^  u    i;3   =  V2iJ>3  .  (2.2) 

The  quasi-geostrophic  vorticity  equation  in  pressure 
coordinates  with  horizontal  eddy  diffusion  is 


ft  v'*  +  'kxv*.v(v2*)  +  bo  |i  -  f0  |a  =  ahv\  .   (2.3) 


Observe  that  when  u  =  0  and  AH  =  0,  (2.3)  becomes  the  baro- 
tropic  vorticity  equation  which  is  sufficient  to  describe 
the  main  instability. 


10 


Apply  this  equation  at  levels  1  and  3  which  gives 


It    V\   +  ""VVVCVV    +    go   53T  "   fo   if  =    AHV^*1      ■    (2^> 


ft   V   *3   +  lkxV^3-V(V  V    U0    9T-fo  -V"  =  AHV  *3 


(2.5) 


If  it  is  desired  to  stop  at  some  intermediate  point,  say  the 
tropopause,  rather  than  to  include  the  entire  atmosphere, 
smaller  layers  could  be  used. 

Next  consider  the  quasi-geostrophic  first  law  of  thermo- 
dynamics in  the  form 


2 

__l  +  DcxV*.(5l)  +  -a,  =  _ —L        ,  (2.6) 

o 


where 


"^     [  H+f  ^   ,  (2.6a) 

P  g  P 


and  T  is  the  temperature  obtained  from  the  standard  atmosphere. 

The  heating  follows  the  formulation  of  Ooyama  (1964), 
Charney  and  Eliassen  (1964),  and  Kuo  (1965),  i.e.,  the 
condensation  heating  is  proportional  to  the  convergence  of 
moisture  in  the  Ekman  layer.   The  form  used  was  given  by 
Ogura  (1964) 
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C  T  — 

Q  =  i  _E_  ii  n  w   ,  (2.7) 


where  n  is  a  non-dimensional  parameter.   This  relationship 
is  used  in  both  rising  and  subsiding  air  since  a  Fourier 
expansion  in  x  is  used;  therefore  the  amplitude  of  heating 
is  divided  by  two. 

When  (2.6)  is  applied  at  level  2  the  result  is 


2 
~                ^-i+^q             ApcfpGOp    -y  KQp 
ft  (^-^3)  +  'kxV(  x2  ^-Vd^-iJ^) 1-=-  =  — 2f-^  .    (2.8) 


Define  the  following  quantities 


y\>1   +    i(>3 


*M  =    2      ■  (2-9) 


ip1  -  ip. 
2~ 


*m  =        O    j     ,  (2.10) 


which  implies  that 


*1  =  %  +  ^T   »  (2,11) 


^3  =  ^M  "  *T   *  (2*12) 


Here  ty^   is  proportional  to  the  layer  thickness  and  is  there- 
fore a  measure  of  the  mean  temperature  in  the  layer.  Use 
these  definitions,  add  (2.4)  and  (2.5)  and  divide  the  result 
by  two,  which  gives 
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|^  V2^M   +  IkxV^-VCV2^)    +  IkxVijy  V(V2^) 
3*M        fowH  k 


For  the  second  forecast  equation,  subtract  (2.5)  from  (2.4) 
and  eliminate  co„  with  the  use  of  (2.8)  which  gives 


|^  (V2-y2)^T  +  lkxV<J,M-V(V2-u2)^T  +  IkxViJyVtV2^) 

+  3o  83T  +  W  =   -2^  +  V\  >  (2'lM> 

o 


where  ~ 

y2  =  — ^—  (2.14a) 

P  a2 


These  are  the  prediction  equations  for  the  model.   To 
linearize  the  equations,  the  flow  is  separated  into  an  east- 
west  current  that  varies  only  in  y,  and  a  small  departure 
from  this  flow.   Thus  it  is  possible  to  treat  waves  in  x 
independently  and  the  fields  can  be  defined  as  follows 

ipM  =  E(y)  +  A(y,t)cos  kx  +  B(y,t)sin  kx   ,         (2.15) 

^T  =  F(y)  +  C(y,t)cos  kx  +  D(y,t)sin  kx   ,         (2.16) 
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where  k  is  the  x  wave  number.   The  coefficients  A  through  D 
are  Fourier  amplitudes  of  the  disturbance  and  E  and  F  are 
the  coefficients  of  the  mean  zonal  fields. 

Substitute  the  expressions  for  ik,  and  ^„  into  (2.13), 
separate  the  various  sine  and  cosine  terms,  and  neglect  all 
products  of  the  quantities  A  through  D.   Equate  the  coeffi- 
cients of  the  cosine  terms  which  gives 


3   ,32A  .,2s      .  r3E  32B  ,  3F  32D    33E  D    93F  ~ 

TTTT      ( O       ~      Ak       )  =  k      [S ST     +  ~ 7T  ~      5"     B  -      5"     D 

at      »    2  ay    ~    2  3y    ,.    2         ^3              <n    3 

9y  J    3y  J   3y         3y             3y 


2 

-   k2(f|  B   +   ||  D)]    -    3oBk   -   K(^-  -   k2)(A-C) 


3y 


4  2 

+   AH(LT  -    2k2    ^-|  +   k4A)       ,  (2.17) 

3y  3y 


The  sine  terms  give 


TC  3y2  .  3y  3y2   3y  17   3y3    ay3 


2 

+  k2#A  +  ^C)]  +  B  Ak  -  K(2-«-  -  k2)(B-D) 
3y      3y   /J     0        g  2 


4         2 
+  A  (<L|  _  2k2  2-|  +  k4B)   .  (2.18) 

3y       3y 
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Repeat  the  procedure  for  (2.14)  for  the  cosine  terms  which 
gives 


3   ,32C    n,2            '2x      .  r3E  32D    3F  32B    3E  n,.  2  .  2, 
■^r   ( — 5-  -  Ck   -  Cy  )   =   k  [ *■  + ~- D(k  +  y  ) 

dr  8y  3y  9y    3y  3y    3y 


-  2S  B(k2  -  y2)  -  2-1  D  -  2-5  B]  -  3  Dk 
3y  9y^    3y^      ° 


2  4         2 

+  (K-S)(2-T  -  k2)(A-C)  +  A  (2-C  _  2k2  i_C  +  k4c)   ^ 

3y  3y       3y 


(2.19) 


The  results  for  the  sine  terms  are 


»   (2!D  .  Dk2  _  Dy2}  =  k  [2E  c(k2  +  y2)  +  3F  A(k2  _  y2) 

dt  3y^  3y  ay 


3E  afc  .  3F  ii*.  +  £f  a  +  3%  c]  +  B  Ck 
3y  3y    3y  3y    3y^    3yJ 


2  4  2      4 

+  (K-S)(2-=-  -  k2)(B-D)  +  AH(^-&  -  2k2  i-£  +  k  D)   , 

3y2  H  ay4      3y2 


where 


(2.20) 


K  -  -2-   (A    ,  (2.21) 

2R0     o 
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and 


s  ~  R  "I — : Pn  (-?~"^   7~"   •  (2.22) 

I e i  3p     o   o 


These  are  the  same  equations  used  by  WSN  except  that  in 
their  equations  F  =  0,  n  =  0,  and  A^  =  0 . 

It  is  necessary  to  derive  the  computational  equation  for 
vertical  motion  beginning  with 


?     3^T      dtyn,  9iKp   KQp 

u.  =  -f—   [  -rrr-   +  u  T-!-  +  v   v-^-  +  7^-   ]  ,       (2.23) 
2    Apo    9t     m  8x     m  3y    2f    ' 


where 


u   = ± tl  (2    2k) 

m      f   3y  \****l 


Vm  =    f-9x-   *  (2^25) 

o 


Substitute  equations  (2.24),  (2.25)  and  (2.16)  into  (2.23), 
collect  terms,  and  again  neglect  all  products  of  A  through  D 
which  gives 


2    /    .   r3C    k  ,3E  n    3FD  ^ 
W2  =  ApT  \C0S  kx  [3t  "  f"  (37  D  "  77   B  )] 


i   j   i   r  3D  .  k  /  8E  0    3F  .  %  -i  .  3F  ,    2 

+  sln  kx  [  at  +  T  (a7  c  "  a7  A)]  +  at  +  ~ 

o  o 


(2.26) 
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Both  the  disturbance  and  the  average  vertical  motion  are 
given  by  this  equation.   The  prediction  equations  and  the 
vertical  motion  equations  are  now  in  the  form  to  be  used  in 
the  model. 
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III.   BOUNDARY  CONDITIONS  AND  FINITE 
DIFFERENCING  SCHEME 


The  finite  difference  scheme  that  was  used  is  illustrated 
below  with  a  sample  variable  M. : 


3y    2H  Uii+1   "i-1^  ^,1; 


92M 


M     1 

?  =  "T  (Mi+l  "  2Mi  +  Mi-l}>        (3*2) 


8y    H 


d3N 


-r  =  — ^r(M    -2M    -M  )-(M  -2M   +M    )  1 . 
dy  dti 

(3.3) 


where  i  is  the  grid  index  and  H  is  the  distance  between  grid 
points . 

Centered  time  differences  were  used  for  all  quantities 
except  those  involving  friction  and  heating.   The  friction 
and  heating  terms  were  computed  at  time  (t  -  At).   Initially 
a  forward  time  step  was  used.   The  second  order  equations 
for  time  tendencies  were  solved  by  the  exact  method  of 
Richtmyer  (1957,  p.  101). 

Yanai  and  Nitta  (1968a)  studied  the  problem  of  using 
finite  difference  approximations  to  solve  dynamic  instabil- 
ity problems  of  non-divergent  barotropic  currents.   They 
showed  that  for  either  symmetric  or  ant i- symmetric  zonal 
currents,  the  exact  boundary  conditions  can  be  replaced  by 
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placing  rigid  boundary  conditions  at  a  distance  equal  to 
one-half  the  width  of  the  shearing  wind  belt.   Therefore 
rigid  boundaries  are  placed  at  y  =  -W/2  and  y  =  W/2,  where 
W  is  the  total  width  of  the  computational  field.   Except 
for  the  E  and  P  fields,  all  variables  have  the  boundary 
condition  A=B=C=D=0aty=±  W/2.   Second  derivatives 
that  are  required  on  the  boundary  are  set  equal  to  zero. 
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Fig.  2  Mean  wind  profile  U  =  -U  tanh  y/y( 
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IV.   WIND  PROFILE  AND  INITIAL  CONDITIONS 

WSN  investigated  the  structure  of  tropical  disturbances 
using  a  similar  model  and  the  wind  profile  (U  =  -U  tanh  y/y  ) . 
It  is  proposed  to  expand  on  that  paper  to  include  a  sinusoidal 
profile  for  an  initial  disturbance  and  a  thermal  wind  field 
which  will  be  a  mirror  image  of  the  mean  wind  field  used 
by  WSN. 

If  the  waves  are  to  grow,  the  zonal  wind  profile  must 
be  barotropically  unstable.   Charney  (1963)  suggested  that 
a  shear  zone  could  be  formed  by  the  transport  of  angular 
momentum  from  different  regions.   This  would  occur  when  the 
Intertropical  Convergence  Zone  is  shifted  away  from  the 
Equator.   The  unstable  profile  chosen  to  represent  this  is 


Ul  =  U3  =  ~Uo  tanh  y/yo'  (4>1) 


which  is  illustrated  in  Figure  2.   To  represent  a  thermal 
wind  field  the  unstable  profile  is 

\J1   =  -U.  =  Uq  tanh  y/yQ.  (4.2) 

Lipps  (1970)  has  used  the  profile  represented  by  (4.1) 
in  an  inviscid  study  of  barotropic  instability  in  the  tropics 
This  wind  field  is  anti-symmetric  with  respect  to  the 
Intertropical  Convergence  Zone,  located  at  y  =  0. 
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For  the  experiments  which  have  the  wind  profile  (4.1), 
a  range  of  non-dimensional  wave  numbers  (ky  )  from  0.10 
to  3^00  will  be  used.   This  wide  range  of  wave  numbers  is 
necessary  because  it  is  uncertain  which  wave  number  will 
be  the  most  unstable  mode.   Those  experiments  with  the  (4.2) 
wind  profile  will  be  run  over  a  range  of  0.10  to  1.20.   One 
experiment  with  the  latter  wind  profile  containing  CISK 
heating  was  run  over  a  wave  number  range  of  0.10  to  3.00. 

ik.  is  introduced  into  the  model  as  a  sinusoidal  distur- 
bance in  y  and  rp      grows  from  zero  as  a  result  of  the  Ekman 
friction  term.   CISK  heating  is  added  which  greatly  increases 
the  growth  rates  of  all  fields,  particularly  the  short  waves. 
Horizontal  friction  is  also  added.   This  represents  a  mixing 
of  the  smaller  scale  waves  and  smooths  out  perturbations 
in  the  vertical  velocity  field.   Thus  at  time  t  =  0 


A  =  sin  3ffi  (4.3) 


B  =  0  (4.4) 

C  =  0  (4.5) 

D  =  0  (4.6) 


E  =  fQ  UQ  yQ  ln[cosh  (y/yQ)]   or   0    (4.7) 


P  =  -f  U  y   ln[cosh  (y/y  )]   or  0    (4.8) 
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V.   COMPUTATIONAL  PROCEDURE 

The  following  constants  were  used  In  all  computations 
except  those  related  to  heating  and  friction  which  were 
zero  if  the  heating  and/or  the  friction  are  neglected: 

p       =   100  centibars, 
o 

y       =   400  or  200  or  100  kilometers, 
o 

W       =   4000  kilometers, 

U       =10  meters  per  second, 

o  *  > 

_5 
f       =   5  x  10   per  second, 
o 

a       =0.4  meters  per  second  squared  per 
centibar  squared, 

8       =   0.0  or  2.29  x  10"   per  meter  per  second, 

Ap      =   50  centibars, 

A       =10  meters  per  second  squared, 

5 
Au      =   1  x  10   meters  per  second  squared, 

a       =   22.5  degrees, 

2  —12 

u       =   2.49  x  10    per  meter  squared, 

H  =  50  or  25  or  12.5  kilometers, 

At  =  1.0  or  0.5  or  0.25  hours, 

n  =  2.0  or  1.0, 

80/8p  =  -0.6  degrees  per  centibar. 

A  forecast  period  of  142  days  was  used  so  that  the  most  un- 
stable mode  would  have  reached  the  maximum  growth  rate.   For 
most  cases,  it  was  found  that  the  growth  rate  at  110  days 
was  the  same  as  the  growth  rate  at  the  longer  period. 
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WSN  introduced  a  random  disturbance  into  the  \|;M  field. 
It  was  found  that  by  using  a  sinusoidal  profile  for  the 
disturbance,  similar  growth  rates  were  obtained  except  in 
the  cases  where  heating  was  added.   The  sinusoidal  profile 
was  used  as  the  disturbance  field  for  all  experiments  with 
variations  of  U,  3,  CISK  heating,  horizontal  and  surface 
friction,  and  the  non-dimensional  wave  number  (ky  ).   Growth 
rates  were  then  computed  by  assuming  the  amplitude  could 
be  written  as 


a1  =  aQ  exp(nUo/yQ)  t.^  (5.1) 


a2  =  aQ  exp(nUQ/yo)  t2,  (5-2) 


where  n  is  the  growth  rate  and  a  is  a  constant.   If  the 

&  o 

ratio  ap/a.  is  formed,  the  growth  rate  can  be  shown  to  be 

yo  ln(a2/ai) 

n  =   U  (tp  -  O   '  (5*3) 

o  d  i 


The  dimensional  doubling  time  is 


yrt  In  2 
o 
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VI.   RESULTS 

A  series  of  experiments  were  conducted  over  a  wide 
variety  of  conditions.   The  total  distance  over  the  grid 
was  ^000  km.   A  grid  mesh  size  of  50  km  and  a  time  step  of 
1.0  hours  was  used  for  all  experiments  except  two  cases. 
In  one  case  a  grid  mesh  of  25  km  with  a  0.5  hour  time  step 
was  used  and  in  the  other  case  a  grid  mesh  of  12.5  km  and 
a  0.25  hour  time  step  was  used. 

The  more  important  experiments  are  given  in  Table  1  as 

a  function  of  U,  $,  heat,  friction,  and  ky  .   The  growth 

rate  n  is  also  given.   The  table  shows  whether  or  not  CISK 

heating  and  horizontal  friction  were  included  for  each 

experiment.   All  experiments  given  in  Table  1  included 

surface  friction.   For  each  experiment,  calculations  were 

made  to  determine  the  non-dimensional  growth  rate  as  a 

function  of  the  non-dimensional  wave  number  ky   so  that  the 

J  o 

results  could  be  compared  to  those  obtained  by  WSN.  The 
value  of  ky  that  was  associated  with  the  maximum  growth 
rate  is  shown  in  the  table. 

It  should  be  mentioned  that  the  value  used  for  the  non- 
dimensional  heating  parameter  x\   was  too  large.   A  more 
realistic  value  for  n  would  have  been  1.0  instead  of  2.0. 
Test  runs  were  made  for  experiments  1  through  4  with 
H  -  1.0  and  in  all  cases,  the  growth  rates  were  reduced  to 
about  one-third  the  value  obtained  with  the  larger  value 
of  n. 
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Table  1.   Compilation  of  numerical  experiments 


Exp 

U 

B 

n 

Heating 

Friction 

ky 
Jo 

1 

°1  =  "3 

0 

0.770 

Yes 

No 

0.40 

2 

u1  =  u3 

2.29X10"11 

0.730 

Yes 

No 

0.40 

3 

Ul  =  U3 

0 

0.520 

Yes 

Yes 

0.40 

4 

°i-*3 

2.29X10"11 

0.520 

Yes 

Yes 

0.40 

5 

Ul  =  ~U3 

0 

0.283 

No 

No 

0.65 

6 

u,  =  -u3 

0 

0.790 

Yes 

No 

1.10 

Experiment  1  was  run  with  the  (4.1)  wind  profile  and 
included  CISK  heating  and  surface  friction.   The  beta  term 
and  horizontal  friction  were  excluded.   Garcia  (1956)  found 
that  the  shear  zone  profile  (4.1)  is  barotropically  un- 
stable when  0  <  ky   <  1.   Betchov  and  Criminale  (1967)  com- 
puted the  eigensolutions  for  this  profile  throughout  this 
range  and  found  the  maximum  instability  at  ky   =0.45  which 
agreed  with  the  results  of  WSN.   The  wave  number  of  maximum 
instability  (0.40)  agrees  quite  closely  with  that  value  but 
the  growth  rate  (n  =  0.770)  is  much  greater  than  the  value 
obtained  by  WSN  (n  =  0.153)  even  when  the  smaller  value  of 
r\   was  used.  Figure  3  contains  the  growth  rate  profile  for 
the  full  range  of  wave  numbers.   The  maximum  growth  rate 
corresponds  to  a  doubling  time  of  0.42  days  and  a  wavelength 
of  6280  km.   The  sharp  irregularities  between  0.10  and  1.00 


26 


arise  because  the  wave  number  increment  in  that  region  is 
0.10  as  compared  to  0.50  between  1.00  and  3.00.   In  the 
non-heating  cases  of  WSN,  there  was  no  growth  for  ky   >  1, 
but  in  experiment  1  there  is  large  growth  which  must  be  due 
in  part  to  the  CISK  process.   The  Figure  shows  that  the 
growth  is  nearly  independent  of  the  scale.   Figure  4  contains 
the  amplitudes  of  ^_  and  oo?  for  this  experiment.   ip,  and  ty- 
are  not  shown  for  clarity  purposes  as  they  have  the  same 
magnitude  as  $„•      ip-,  is  slightly  larger  and  ip~    slightly 
less  than  iprp  but  the  shape  is  identical.   The  units  are 
arbitrary  since  the  basic  equations  are  linear.   However, 
ojp  is  not  plotted  on  the  same  scale  as  ^„.   All  quantities 
show  double  maxima  with  peaks  at  y  =  ±1.2.   These  results 
are  somewhat  similar  to  those  obtained  by  WSN  without 
heating  but  their  results  showed  the  minimum  between  peaks 
as  only  slightly  less  than  the  maximum  whereas  in  Figure  4, 
the  minimum  is  less  than  one-half  of  the  maximum.   Also  the 
iJj,  ,  4>^t    and  ip^  fields  are  very  close  together  but  WSN 
showed  ip,  to  be  predominate.   The  closeness  of  the  ip  fields 
roughly  implies  that  the  vertical  mean  streamfunction  is 
nearly  zero  so  the  fields  are  basically  of  a  thermal  charac- 
ter.  Figure  5  contains  the  phase  angles  for  iK  ,  ip_,  ^T, 
and  co_  for  this  experiment  as  functions  of  y.   ij>,  and  $„ 
are  very  nearly  coincident  at  the  outer  edges  and  ty-,    =  ^T 
between  y  =  ±2  although  only  i>1   is  shown  in  Figure  5.   iJn 
is  approximately  180  degrees  out-of-phase  from  iK  and  ^T 
which  indicates  the  pressure  trough  at  the  surface  is  replaced 
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Fig.  3   Growth  rates  for  experiment  1. 
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Fig.  4   Amplitudes  of  i^T,  and  u>2  for  experiment  1 
^  and  Up  are  on  different  scales. 
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Fig.  5   Phases  of  i> 
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\p    s    4>T,  and  w2  for 


experiment  1  in  degrees.   Each  dependent 
variable  is  expressed  in  the  form 
Pcos(kx-<5),  where  P  is  the  amplitude 
and  6  is  the  phase  angle. 
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by  a  ridge  at  the  upper  level.   This  is  consistent  because 
the  Ekman  effect  and  heating  are  tied  to  the  field  at  level 
3  which  produces  heating  above  surface  low  pressure  systems. 
Unlike  the  results  found  by  WSN  for  the  adiabatic  case, 
Figure  5  shows  the  horizontal  tilt  of  the  disturbance  to 
be  a  minimum  at  y  =  0  and  a  maximum  at  y  =  ±(1  to  4).   The 
general  sense  of  the  tilt  is  northeast  but  larger  than  that 
observed  in  easterly  waves.   Also  the  tilt  of  the  cjo?  field 
fluctuates  considerably.   The  small  perturbations,  particu- 
larly in  the  Wp  field,  suggest  that  there  is  not  a  clearly 
defined  unstable  mode  in  y.   This  does  not  agree  with  the 
results  of  WSN  which  had  no  heating.   It  is  interesting 
to  note  that  the  fields  are  very  smooth  but  have  large 
gradients  except  for  u)p  which  shows  a  radical  behavior. 
Experiment  2  is  the  same  as  experiment  1  except  the 
beta  term  is  included.   The  results  in  general  are  the  same. 
Figure  6  shows  the  growth  rate  profile.   The  magnitude  of 
the  growth  rates  is  almost  the  same  although  a  small  decrease 
in  the  maximum  is  noted  (from  0.770  to  0.730)  and  the 
perturbations  betwen  ky   =  0.10  and  1.00  have  been  somewhat 
smoothed.   The  maximum  growth  rate  corresponds  to  a  doubling 
time  of  0.44  days  and  a  wave  length  of  6280  km.   Figure  7 
contains  the  amplitudes  for  experiment  2  as  a  function  of 
y.   As  in  experiment  1,  only  ipm  and  Up  are  shown  with  i>^ 
and  to  on  different  scales.   Again  there  is  a  double  maxima 
at  y  =  ±1.2  but  the  profile  is  asymmetric  with  the  larger 
maximum  in  the  shear  zone  to  the  north.   This  asymmetry  is 
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Fig.  6   Growth  rates  for  experiment  2. 
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Fig.  7   Amplitudes  of  ^T.and  co?  for  experiment  2 
^T  and  o)p  are  on  different  scales. 
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Fig.    8      Phases    of   iK  ,    i|>~,    ^T>    and   w2    for 

experiment  2  in  degrees.   Each  dependent 
variable  is  expressed  in  the  form 
Pcos(kx-6),  where  P  is  the  amplitude 
and  6  is  the  phase  angle. 
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expected  since  the  absolute  vorticity  gradient  is  not  an 
odd  function  of  y.   WSN  found  that  the  shift  of  symmetry  was 
to  the  south  while  Figure  7  shows  a  shift  to  the  north. 
The  phase  relationships  of  ip,  ,  ^_,  ip_ ,  and  co_  shown  in 
Figure  8  are  very  similar  to  those  in  experiment  1  and 
again  only  ^»  is  shown  between  y  =  ±2. 

Experiment  3  was  run  with  the  ( ^ .  1 )  wind  profile  and 
included  CISK  heating,  surface  friction,  and  horizontal 
friction.   The  beta  term  was  excluded.   The  horizontal 
friction  was  added  in  an  attempt  to  smooth  out  the  vertical 
velocity  (oj?)  pattern.   Figure  9  shows  the  growth  rate 
profile  for  experiment  3.   The  maximum  growth  rate  has  been 
reduced  considerably  (from  0.770  to  0.520)  and  the  growth 
falls  off  at  ky   >_  2.00.   The  maximum  growth  rate  corresponds 
to  a  doubling  time  of  0.62  days  and  a  wavelength  of  6280  km. 
Note  that  there  are  no  perturbations  in  the  profile.   Figure 
10  contains  the  amplitudes  of  ip™  and  Wp  for  this  experiment. 
As  in  experiment  1,  only  4>T  and  Wp  are  shown  and  are  on 
different  scales.  'if>,  is  slightly  larger  and  i/>~  slightly 
smaller  than  ty™.   It  is  interesting  to  note  that  the  largest 
amplitudes  are  located  near  the  boundaries  for  all  fields 
whereas  in  the  case  without  friction  (experiment  1),  the 
maximum  amplitudes  are  in  the  central  regions.   Perhaps  the 
large  scale  eddies  favored  by  the  horizontal  friction  are 
heavily  damped  by  the  horizontal  shear  between  y  =  ±1. 
Figure  11  contains  the  phases  of  ij>,  ,  ty~   and  o)p.   ip_  is  not 
shown  but  it  is  everywhere  equal  to  if>,  .   Again  ^T  and  i>~ 
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Fig.  9   Growth  rates  for  experiment  3- 
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Fig.  10   Amplitudes  of  ^T  and  oj2  for  experiment  3 
i^m  and  Wp  are  on  different  scales. 
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Fig.  11   Phases  of  iK  ,  i|>-,  and  o)p  for  experiment 

3  in  degrees.   Each  dependent  variable 
is  expressed  in  the  form  Pcos(kx-6), 
where  P  is  the  amplitude  and  6  is  the 
phase  angle. 
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are  out-of-phase  but  all  fields  are  smoother  than  in  experi- 
ment 1.   Vfhere  there  are  peaks  in  the  amplitude  fields,  the 
phase  angles  of  u)_  and  i^_  are  the  same  which  tends  to 
suggest  that  warm  air  is  rising  and  cold  air  is  sinking. 
Again  there  is  an  eastward  tilt  in  the  i>   fields.   The  wp 
field  tilt  is  opposite  to  the  ip  fields  but  it  does  correlate 
with  the  ty-    field  where  the  amplitude  is  a  maximum.   It 
appears  the  scale  in  y  is  about  2y  . 

Experiment  4  is  the  same  as  experiment  3  but  the  beta 
term  is  included.   The  beta  term  has  very  little  effect 
except  to  make  the  solutions  asymmetric.   Figure  12  shows 
the  growth  rate  profile  which  is  the  same  as  Figure  9 
except  for  the  slight  perturbation  between  the  wave  numbers 
0.40  and  0.90.   The  maximum  growth  rate  corresponds  to  a 
doubling  time  of  0.62  days  and  a  wavelength  of  6280  km. 
Figure  13  contains  the  amplitudes  of  ^  and  oj?  .   As  in 
experiment  3S    $-,    and  ip«  are  not  shown.   The  only  change  noted 
is  a  slight  increase  in  the  maximum  to  the  north  and  a 
decrease  in  the  maximum  to  the  south.   Figure  lh   contains 
the  phases  of  ip,  ,  i|>_  and  (jo„  and  again  ip-  is  everywhere  equal 
to  ip,  .   The  oj_  and  \\>^   correlation  is  the  same  as  in  experiment 
3. 

Experiment  5  was  run  with  (4.2)  wind  profile  and  surface 
friction.   Heating  and  horizontal  friction  were  excluded. 
This  profile  is  a  mirror  reflection  of  the  U,  profile,  that 
is  -U,.   This  profile  is  similar  in  general  to  the  wind 
profile  in  the  Intertropical  Convergence  Zone  where  there 
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is  cyclonic  shear  in  the  lower  region  and  anticyclonic 
shear  in  the  upper  region.   Since  the  ITCZ  is  driven  by 
heating  and  convergence,  this  is  the  type  of  mean  wind 
pattern  that  is  expected.   Three  cases  were  run  although  only 
one  is  shown  in  the  table.   Case  1,  which  is  summarized  in 
Table  1,  had  a  grid  mesh  of  50  km,  a  yQ  value  of  400  km  and 
a  time  step  of  1.0  hours.   Case  2  used  a  grid  mesh  of  25  km, 
a  y   value  of  200  km  and  a  time  step  of  0.50  hours.   Case  3 
used  a  grid  mesh  of  12.5  km,  a  yQ  value  of  100  km  and  a 
time  step  of  0.25  hours.   Case  1  was  run  over  the  range  of 
wave  numbers  from  0.10  to  1.20  and  the  maximum  unstable 
mode  was  found  at  kyQ  =  0.65.   Cases  2  and  3  were  run  only 
at  wave  number  O.65.   Figure  15  shows  the  growth  rate  profile 
for  experiment  5,  case  1.   The  maximum  growth  rate  at  kyQ  =0.65 
was  0.283.   The  doubling  time  was  1.13  days  with  a  wavelength 
of  3465  km.   For  case  2,  the  growth  rate  obtained  was  0.188 
with  a  doubling  time  of  0.85  days  and  a  1930  km  wavelength. 
For  case  3  a  growth  rate  of  0.162  was  found  with  a  doubling 
time  of  0.50  days  and  a  967  km  wavelength.   For  cases  2  and 
3,  ky   =  0.65  may  not  be  the  maximum  growth  rate.   This  shows 
that  as  the  scale  (yQ)  is  decreased,  the  baroclinic  process 
is  suppressed.   This  is  consistent  with  baroclinic  instab- 
bility  theory.   It  is  noted  that  in  case  3,  the  growth  rate 
is  the  same  as  that  obtained  by  WSN  in  their  experiment  1. 
Figure  16  contains  the  amplitudes  of  lj^,  ^  and  u>2  for  case 
1.   For  clarity,  i|>-  is  not  shown  but  it  is  slightly  less 
than  «J>, .   These  amplitudes  are  very  similar  to  those  of  WSN 
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Pig.  12   Growth  rates  for  experiment  4 
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Fig.  13   Amplitudes  of  ^T  and  wp  for  experiment 
4.  \p~   and  Wp  are  on  different  scales. 
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Fig.  Ik      Phases  of  i> 
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\p    ,  and  a)p  for  experiment 


h   in  degrees.   Each  dependent  variable 
is  expressed  in  the  form  Pcos(kx-5), 
where  P  is  the  amplitude  and  6  is  the 
phase  angle. 
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in  their  Figure  4 A.   The  amplitudes  for  cases  2  and  3  are 
not  shown  but  they  are  much  the  same  as  Figure  16.   Figure 
17  contains  the  phases  of  vJj.  ,  i|^,  \Jv,,  and  w   as  functions 
of  y.  ^m  and  w   are  approximately  165  degrees  out-of-phase . 
This  implies  a  baroclinic  conversion  of  kinetic  energy. 
However,  the  tilt  is  opposite  to  the  shear  which  indicates 
a  barotropic  conversion  of  kinetic  energy.   The  tilt  is 
westward  with  height  north  of  the  ITCZ  and  eastward  with 
height  south  of  the  ITCZ  which  is  consistent  with  rising 
motion  and  unsettled  weather  that  is  often  observed  east 
of  the  surface  trough  axis  in  easterly  waves  north  of  the 
ITCZ.   These  results  are  also  consistent  with  those  of 
WSN. 

Experiment  6  was  run  with  the  (4.2)  wind  profile,  CISK 
heating  and  surface  friction.   The  beta  term  and  horizontal 
friction  were  excluded.   Figure  18  shows  the  growth  rates 
for  this  experiment  with  a  maximum  unstable  mode  at  ky   =1.10 
corresponding  to  a  doubling  time  of  0.41  days  and  a  wave- 
length of  2285  km.   The  growth  rate  of  0.790  is  similar  to 
that  of  experiment  1  and  again  the  large  value  was  due  to 
the  large  value  of  n.   The  perturbations  between  wave  numbers 
0.10  and  1.10  are  even  larger  than  in  experiment  1  indicating 
that  there  is  no  clearly  defined  unstable  mode  in  y.   The 
amplitudes  have  a  pattern  similar  to  Figure  h.      The  phases 
for  this  experiment  tilted  northeast  as  in  the  other  cases. 


W 


.00  -i 


.70- 


.60- 


.50- 


n  .40 


.30- 


.20  - 


J0- 


Fig.  15   Growth  rates  for  experiment  5,  case  1 
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10 


Fig.  16   Amplitudes  of  ty    ,  ^_,  and  to?  for 

experiment  5,  case  1.   The  quantities 
<K  and  ^„  have  the  same  scale;  u>~   is 

on  a  different  scale. 
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5  n 


»  0 


Pig.  17   Phases  of  iK,  ty~,    ipT,  and  w2  for 

experiment  5,  case  1  in  degrees.   Each 
dependent  variable  is  expressed  in  the 
form  Pcos(kx-6),  where  P  is  the 
amplitude  and  6  is  the  phase  angle. 
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Fig.  18   Growth  rates  for  experiment  6 
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VII.   CONCLUSIONS 

In  this  thesis,  numerically  computed  bare-tropically  un- 
stable disturbances  are  compared  with  tropospheric  waves 
that  occur  in  the  tropical  easterlies.   A  two-level  quasi- 
geostrophic  model  is  used  with  an  attached  Ekman  layer.   The 
most  unstable  solutions  are  obtained  numerically  by  the 
initial  value  procedure.   Two  basic  shear  zone  profiles  are 

used:  (1)  U.  =  U  = -U   tanh  y/y   and  (2)  U_  =  -U0  =  U  tanh  y/y  . 
1   3    o      J    J o         1    3   o      JJo 

Calculations  were  made  in  which  CISK  heating  and  horizontal 
friction  were  added. 

Certain  features  of  the  predicted  wave  structures  do 
correspond  to  observed  tropical  wave  disturbances.   The  tilt 
of  the  disturbance  in  a  sense  opposite  to  the  horizontal 
wind  shear  that  was  observed  by  Wallace  and  Chang  (1969)  ^s 
present  in  all  the  experiments.   All  experiments  which  contain 
heating  have  the  warm  structure  which  was  observed  by  Chang 
et' al.  (1970).   The  presence  of  rising  motion  to  the  east 
of  the  wave  axis  which  was  noted  by  Yanai  and  Nitta  (1967) 
is  observed  in  experiment  5. 

The  growth  rates  obtained  with  heating  are  unrealistically 
large,  but  it  was  shown  that  these  could  be  greatly  reduced 
by  choosing  a  smaller  value  of  the  parameter  n .   It  was  also 
noted  that  when  heating  is  present  there  is  no  truely 
preferred  scale  of  instability.   This  has  been  observed  by 
Chang  (1971).   It  was  pointed  out  in  the  Introduction  that, 
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while  mature  tropical  wave  disturbances  may  be  driven  by 
CISK  heating,  they  may  originate  from  barotropic  instability 
(Bates  1970).   If  so,  then  the  disturbance  should  retain 
certain  features  (latitude  structure,  wavelength,  etc.) 
of  the  barotropically  unstable  modes.   This  would  provide 
a  selection  process  which  does  not  occur  if  CISK  heating 
is  present  from  the  start. 

The  studies  of  this  thesis  should  be  continued  but  a 
more  complicated  model  is  needed.   Proper  parameterization 
of  the  heating  term  is  necessary.   A  primitive  equation 
model  which  includes  a  better  treatment  of  the  boundary 
layer  should  be  used. 
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